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We explore the statistical properties of the Casimir-Polder potential between a dielectric sphere 
and a three-dimensional heterogeneous medium, by means of extensive numerical simulations based 
on the scattering theory of Casimir forces. The simulations allow us to confirm recent predictions for 
the mean and standard deviation of the Casimir potential, and give us access to its full distribution 
function in the limit of a dilute distribution of heterogeneities. These predictions are compared 
with a simple statistical model based on a pairwise summation of the individual contributions of 
the constituting elements of the medium. 
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I. INTRODUCTION 

Materials placed in a close vicinity to each other mod¬ 
ify the modes of the electromagnetic field. This results in 
a change of the vacuum energy, which eventually mani¬ 
fests itself as a net force known as the Casimir force Hi- 
The Casimir force has been the subject of a number of 
experimental investigations at object separations ranging 
from tens of nanometers to a few micrometers. Starting 
with the experiments by Lamoreaux jjjj and Mohideen Q , 
the Casimir effect has experienced an enormous increase 
in experimental activities in recent years H-®. 

Theoretical approaches to the Casimir force are usually 
built on an effective medium description of the interact¬ 
ing materials. Within such an approximation, the lo¬ 
cal details of the materials’ microstructure are neglected 
and the objects are described by macroscopic, spatially- 
independent dielectric constants. While the effective 
medium description is in general quite satisfactory for 
describing dense materials that indeed look homogenous 
at the typical scales of the Casimir force, this is not neces¬ 
sarily the case for strongly heterogeneous (“disordered”) 
media that are made of many constituting elements 
(“scatterers”) well separated from one another. Exam¬ 
ples of such heterogeneous systems include nanoporous 
materials [21], clouds of cold atoms [2^1 an d, in a slightly 
different context, corrugated surfaces |23lI24{. 

From a theoretical viewpoint, interaction phenomena 
involving strongly heterogeneous materials have been lit¬ 
tle studied. Seminal works on that subject considered the 
thermal Casimir interaction between slabs made of piled 
layers separated from random distances (one-dimensional 
disorder) ®, ®|. The question of disorder was also ad¬ 
dressed recently 27] in the context of the Casimir-Polder 
(CP) interaction [28j] between a sphere and a plate ®]. 
In a recent work finally, the CP interaction between a 
dielectric sphere (or an atom) and a three-dimensional 
disordered dielectric material was also investigated f30j . 
This is the scenario we consider in the present paper. 

When a probe sphere or an atom interacts with a spa¬ 
tially heterogeneous material such as a semi-infinite dis¬ 
ordered medium, the CP potential naturally fluctuates in 


space. In other words, the Casimir interaction depends 
on the specific statistical realization of the disorder. A 
shared conclusion of Refs. [2514271 . ® is that when the 
two objects are far enough from each other, precisely 
when the distance between them is large compared to 
the typical separation between two heterogeneities, the 
value of the Casimir potential from a realization to an¬ 
other is well captured by its configuration average, which 
coincides with the prediction of the effective medium pre¬ 
scription. In strong contrast, at smaller distances fluctua¬ 
tions of the potential become larger than its mean, which 
is consequently no longer representative. In practice, this 
conclusion is crucial for measurements of quantum reflec¬ 
tion [2ll l3ll - [34| , and more generally for any measurement 
of the Casimir force involving heterogeneous materials. 

In our previous work [30), we developed an exact math¬ 
ematical treatment of the fluctuations of the CP interac¬ 
tion between a dielectric sphere and a dilute disordered 
dielectric medium, and applied it to the calculation of 
the mean value of the CP potential and of its standard 
deviation. In this paper, we consider the same geome¬ 
try (recalled in Sec. Et, for which we perform extensive 
numerical simulations of the CP potential. The results 
of these simulations confirm the predictions of [3(J (Sec. 
ED, and additionally allow us to compute the full proba¬ 
bility distribution of the CP potential which, for a given 
distribution of the scatterers, does not depend on the 
microscopic properties of the latter. In a second time 
(Sec. T71 . we present a simple statistical model based 
on a pairwise summation of the individual contributions 
of the scatterers, and confront it with the simulations. 
Concluding remarks are collected in Sec. |V] 

II. MEAN AND STANDARD DEVIATION OF 
THE CASIMIR-POLDER POTENTIAL 

We address the CP interaction between a probe di¬ 
electric sphere (placed in vacuum) of static polarizabil¬ 
ity ao (here and in the following, polarizabilities are ex¬ 
pressed in SI units divided by eo) and a semi-infinite, 
three-dimensional disordered medium consisting of a col¬ 
lection of many scatterers, as illustrated in Fig. [T] We 


2 


denote by z the distance between the sphere and the sur¬ 
face of the disordered medium. For definiteness, in this 
paper we restrict our discussion to the retarded regime of 
the Casimir interaction where z much exceeds the reso¬ 
nance wavelength Ao of the probe sphere (the treatment 
of the opposite limit z -C Ao is analogous). Scatterers 
are also modeled by dielectric spheres of size a and of 
static polarizability a s . Throughout the paper, we as¬ 
sume that they are uniformly distributed in space with 
density n, and we consider the case of a dilute disor¬ 
dered medium, such that the average distance between 
the scattering spheres is larger than their size, na 3 -C 1 . 
This is the typical regime where the statistical fluctu¬ 
ations of the CP potential are the largest ^3]. In the 
opposite limit na 3 ~ 1 of scatterers very close to each 
other, the approach developed below does not apply but 
we expect the statistics of the Casimir-Polder potential 
to be approximately captured by a model where the atom 
interacts with a rough surface [351 ] . 





CP potential and its mean, was found to be (for z Ao): 


\l&U*{z) 
\U(z) | 


a i 
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with a\ ~ 0.7. Equation fljZJ indicates that U(z), i.e. the 
prediction of the effective medium theory, is well rep¬ 
resentative of U{z) only when z n -1 / 3 . At smaller 
scales, 7 becomes larger than unity and U(z) no longer 
provides a trustful estimation of the interaction. 


III. NUMERICAL SIMULATIONS 
A. Methodology 

We now propose to investigate the statistical properties 
of the CP potential from exact numerical simulations in 
the geometry of Fig. [I] For this purpose, we proceed as 
follows. We generate an ensemble of N dielectric spheres 
of radius a and frequency-dependent permittivity e(w), 
uniformly distributed in a cube of side L. This system 
constitutes a disordered medium of average density n = 
NL~ 3 . An additional probe sphere is placed above this 
cube, at a distance z to the center of one face, as in Fig. |T| 
Denoting by U < ' N+1 ' > (z) the total, internal Casimir energy 
between the N - 1-1 spheres [if|, the CP interaction U(z) 
is by definition 

U(z) = U {N+1 \z) - U {N+1 \z -A oo) 

= £/ (iV+ 1 ) (z)-t/ (JV) . (3) 


FIG. 1. (color online) We consider the Casimir-Polder in¬ 
teraction between a dielectric sphere (placed in vacuum) and 
a semi-infinite disordered medium. The disordered medium 
consists of a collection of dielectric spheres (size a, density n) 
whose positions are uniformly distributed in space. 

In [30| , the question of fluctuations in the limit na 3 <C 
1 was tackled with the help of a statistical description 
of the disordered material, in which the CP potential 
U(z) becomes a random variable. Its mean, U(z), and 
its variance, 5U 2 {z ), were calculated from an exact treat¬ 
ment of radiation-matter interaction, based on the scat¬ 
tering approach to Casimir forces [36J, |37| combined with 
a diagrammatic description of radiation scattering off the 
disordered medium [38|, |39| . In the limit z Ao, the fol¬ 
lowing expression for the mean was found [30| : 
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U(z) = - na s U*(z ), (1) 

where U*(z ) = — 3ao^c/(327r 2 ^ 4 ) is the Casimir poten¬ 
tial between the probe sphere and a perfect mirror. As 
announced, the result ([!]) coincides with the prediction of 
an effective medium description where the probe sphere 
interacts with an homogeneous surface of relative per¬ 
mittivity e = 1 + na s . The amplitude of fluctuations, 
quantified by the ratio 7 of the standard deviation of the 


The strategy thus consists in calculating the interaction 
energy as the difference between the internal energies of 
N + 1 and N spheres. Within the scattering formal¬ 
ism [36;, '37} , the total Casimir energy between N spheres 
is given by I ID 

h f°° 

UW = — dwlogdet (MM- 1 ) . (4) 
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M is a block-square matrix of dimension N with the fol¬ 
lowing structure: 


M = 


( 1 \ 

Ri Ti _> 2 Ti^jv 

T2-J.I R-2 1 ■ ’ ' T 2 _ > at 


( 5 ) 
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Tjv_>2 


R 


-1 

’N 




The diagonal blocks of M are the inverse of the spheres’ 
reflection operators R,. The (i, j) off-diagonal block of M 
contains the translation operator Tj_>y, which relates an 
outgoing spherical wave centered on 17 to an incoming 
spherical wave centered on 17 |4C|. Finally, M ^, 1 is the 
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(color online) Absolute value of the mean CP po¬ 
tential between a dielectric sphere and a semi-infinite disor¬ 
dered medium, as a function z. Dots are the results of exact 
numerical simulations and the dashed line is the theoretical 
prediction ©• 
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block-diagonal matrix diag(Ri,... Rat). For the simula¬ 
tions, we express the scattering and translation operators 
Ri and T,;_yj in a basis of spherical vector waves | £mP), 
with £ > 1, —£ < m < £ and P = {E,M}. In this basis, 
the matrix elements of R, are given by the standard Mie 
scattering amplitudes 0] . We compute these amplitudes 
without any approximation, taking into account the full 
multipole expansion. Finally, we evaluate the matrix el¬ 
ements of using the formalism of Ref. 0] • 



FIG. 3. (color online) Relative fluctuations of the CP po¬ 
tential as a function of nz 3 . Dots are the results of exact 
numerical simulations and the dashed line is the theoretical 
prediction ©• 
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B. Mean and standard deviation 

Making use of the approach described in Sec. IIII A1 
we compute the CP potential U(z) between a dielectric 
sphere of radius a = 10 nm and a disordered medium con¬ 
sisting oi N = 32 other, identical spheres with the same 
radius a, uniformly distributed in a cube of side L = 
12 n m. For these parameters, the disordered medium is 
dilute, na 3 -C 1, and we are effectively describing the ge¬ 
ometry of a semi-infinite bulk system as long as z <C L. 
For definiteness, we give to all the spheres the permittiv¬ 
ity of silicon: e(w) = 1 — [e(0) — 1]wq/(u; 2 — Wq + 
where wo = 27tc/Ao with Ao = 295nm, e(0) = 11.6 and 
7 = 0.03u;o 0. Figure [2] displays the absolute value 
of the disorder-averaged CP potential, C/(z), computed 
with these parameters, for several values of 2 (red dots). 
Each dot is obtained by generating thousands (typically 
between 3000 and 8000 depending on z) of disorder re¬ 
alizations, computing U(z) for each of them and finally 
averaging the results. In Fig. [2] we also show the theo¬ 
retical prediction (JT]) (dashed line), which is in very good 
agreement with the numerics. The small disagreement 
visible at large z stems from deviations to the geometry 
of the semi-infinite medium: when z becomes of the or¬ 
der of L/2 , the probe sphere starts to be sensitive to the 
boundaries of the system, and a cross-over toward the 
sphere-cube geometry is expected. We also show in Fig. 
[3] the standard deviation of the CP potential relative to 


FIG. 4. (color online) Scatter plot of the CP potential nor¬ 
malized to its mean value, s = U/U, computed for 5500 dis¬ 
order realizations, for nz 3 = 0.5. 

its mean, 7 , as a function of nz 3 . Red dots are the nu¬ 
merical results, and the dashed curve is the theoretical 
prediction (|2|) . The agreement between theory and nu¬ 
merics is very good, up to small finite-size effects at large 
z. 


C. Probability distribution function 

As was pointed out in [13, the fluctuations of U(z) be¬ 
come significant at distances z < n -1 / 3 , when 7 becomes 
larger than unity, see Eq. © and Fig. [3] This suggests 
that at small distances, the mean U(z) is no longer rep¬ 
resentative of U(z). In order to confirm this picture, we 
compute the full probability distribution function p(s ) of 
the CP potential normalized to its mean, s = U/U , by 
constructing histograms of the numerical data. In Fig. 
© we show as an example a scatter plot of the data ob¬ 
tained for 5500 disorder realizations, for nz 3 = 0.5. The 
associated histogram p(s) is displayed in Fig. © together 
with the histograms corresponding to three other values 
of the parameter nz 3 . A quick look at the distribu¬ 
tions in Fig. [5] forthwith confirms the property already 
outlined by the analysis of 7 : as the sphere gets closer 
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FIG. 5. (color online) Probability distribution function 
p(s = U/U), for four increasing values of nz 3 . Histograms 
(vertical lines) are the results of exact numerical simulations, 
and solid curves are the theoretical prediction, Eq. (1151) . The 
lower-right panel also displays as a dashed curve the Gaussian 
distribution expected in the limit of very large distances, Eq. 

un¬ 


to the disordered medium, the distribution function be¬ 
comes more and more peaked around a value s -C 1, 
corresponding to a CP potential much smaller than its 
mean. In other words, U(z) is no longer a self-averaging 
quantity. Only when nz 3 > 1 does the maximum of the 
distribution approaches s = 1. Such a phenomenon was 
previously observed in the context of the interaction be¬ 
tween plates with one-dimensional disorder miH]. We 
see here that it is a quite general property, not restricted 
to one-dimensional systems. 


D. Sensitivity to microscopic parameters 

To conclude our numerical study, we address the ques¬ 
tion of the sensitivity of p(s) with respect to changes 
in the microscopic properties of the scatterers. For this 
purpose, we perform additional numerical simulations in¬ 
volving scatterers with a different radius a = 1 nm and 
made of a different material with frequency-independent 
permittivity e = 10. We also set the radius of the probe 
sphere to a = 1 nm, but keep the same value of the per¬ 
mittivity of silicon as in the previous section. The distri¬ 
bution p(s ) obtained for these new parameters is shown 
in the right panel of Fig. [6l for nz 3 = 1. For compari¬ 
son, we also display the distribution p(s) computed with 
the parameters of the previous section. No visible differ¬ 
ence is seen between the two histograms, which indicates 
that in the dilute limit p(s) is in fact a function of the 
parameter nz 3 only. In particular, the parameter na 3 is 
irrelevant. This could have been anticipated since in the 
limit na 3 -C 1 of independent scatterers, na 3 enters both 



FIG. 6. (color online) Numerical probability distribution p(s) 
for two different sets of parameters of the scattering spheres, 
for nz 3 = 1 (histograms, vertical lines). No visible difference 
is seen between the two histograms. This is confirmed by the 
theoretical prediction (solid curve), Eq. (1151) . which depends 
on the single parameter nz 3 . 


U and U within the same prefactor, which thus cancels 
out when considering the ratio s = U/U (see Sec. IIVI for 
a general proof). This property is in particular fulfilled 

by the second moment of the distribution, y 2 = SU 2 /U 2 , 
see Eq. ([2]). 


IV. SIMPLE MODEL 

We now develop a simplified statistical description of 
the CP interaction between a dielectric sphere and a 
disordered bulk medium, based on a pairwise summa¬ 
tion (PWS) approximation [Hj]. This approximation de¬ 
scribes the total CP interaction U ( z ) as a sum of the pair 
interaction £ between the probe sphere and each of the 
N scatterers. It has to be distinguished from the pertur¬ 
bative expansion, as it can in principle be used for non 
perturbative pair interactions. In the problem studied in 
this paper however, the validity of the two approxima¬ 
tions is a consequence of the same assumption of a dilute 
disorder ( na 3 <1). 

A. Pairwise summation 

As in the numerical simulations (Sec. Mil) , we consider 
a situation where the distance z much exceeds the sphere 
resonance wavelength Ao (In the opposite limit x < Ao, 
the reasoning follows exactly the same lines). Con¬ 
sequently, the interaction potential between the probe 
sphere and a scatterer located at distance r takes the sim¬ 
ple form £ = —U-j/r 7 , where T 7 is a constant character¬ 
istic of the microscopic properties of the two interacting 
objects. Let us consider a small spherical cap of volume 
dV containing dN scatterers, as illustrated in the left 
panel of Fig. [7] The elementary CP interaction between 
them and the probe sphere is dU = £dN = —TjdN/r 7 . 
As in the simulations, we assume the positions of the 
scatterers to follow a uniform distribution, so the random 
variable dN is Poisson distributed. The m th cumulant of 
































two atoms [45l | : 
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FIG. 7. (color online) Left: parametrization of the geometry 
for the statistical approach of Sec. IIV Al Right: typical (rare) 
disorder realization contributing to the Lifshitz tail: the atom 
is above a large region free of scatterers, of size R^> z. 


dN , K m (dN ), thus fulfills: 

K m (dN) = Ki(dN) = ndV, , (6) 

where n = N/L 3 is the average density of scatterers. 
Within the PWS approximation, the cumulant of the el¬ 
ementary CP potential dU involving dN scatterers reads: 

K m (dU) = K m (dN) 

ndV = dK m (U), (7) 

where we have used the property K m (aX) = a m K m (X) 
in the first equality, and Eq. m in the second. The 
cumulant of the total CP potential between the probe 
sphere and the N scatterers is finally obtained by using 
the parametrization dV = 2nr(r — z)dr , see Fig. [71 and 
integrating over r from r = z to oo. This yields 

(—l) m 2irnT™ 

’ (7m — 3)(7m — 2)z 7m ~ 3 ’ l j 

from which various statistical properties can be deduced, 
as we now discuss. 



B. Mean and standard deviation 


As a preliminary application of our statistical ap¬ 
proach, we propose to re-derive the predictions © and 
©, previously obtained from an exact treatment of the 
radiation-matter interaction. For simplicity we assume 
the scatterers to be spheres of static polarizability op, 
identical to the probe sphere. The coefficient T7 then de¬ 
scribes the large-distance interaction between two iden¬ 
tical spheres. It can be readily evaluated, for instance 
from the Casimir-Polder law [I] for the energy between 

I 


fir 3 r 

pis) = lo*°l 


exp • 


f 2irnz 3 
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To'” t “ 1 


r 2Zhccf ) 
7= (4t r) 3 ‘ 


(9) 


The mean U(z) is by definition the first-order cumu¬ 
lant: U(z) = Ki(z) = —2TrnYi/(20z A )i where we have 
used Eq. © for m = 1. Combining this expression with 
Eq. ©, we recover Eq. ©. The variance SU 2 (z ) is, 
on the other hand, given by the second-order cumulant, 
K 2 (z) = 27mr 2 /(1322: 11 ). Taking the ratio with K\(z), 
we find 


\Jsu 2 (z) ^ i 

\U(z)\ ~ V 33 tt ^3 


( 10 ) 


which is nothing but Eq. ©, with an analytic expression 
of the coefficient a/50/(337t) ~ 0.7. 


C. Probability distribution function 


Let us now derive the probability distribution function 
p(s = U/U) that has been studied numerically in Sec. 
IIII Cl p(s) is given by the inverse Laplace transform 

1 p5-\-ioo 

p(s) = — / e st e^dt, (11) 

Z7TZ Jfi—ioo 

where 5 is greater than the real part of all singularities of 
e v< d\ ip(t ) is the cumulant generating function of s, and 
can be expressed as a power series of the cumulants ©: 


<P(t) = 

m—0 


(■ ~t) m K m (U) 

to! if 71 


( 12 ) 


Making use of Eq. © and of the relation U(z) = 
—27mT7/(202: 4 ) obtained above, we find: 


p(t) 


2irnz 3 

6 


— 1 + e~ T 


2t 3/7 74/ 7 (t) + 3r 2 / 7 7 5/7 (r) 


(13) 

where we have introduced r = 20t/(2irnz 3 ) and where 
lf T ) = fo x q 1 e X dx is the lower incomplete gamma 
function. Since has no singularities in the complex 
plane, we can set i5 = 0 in Eq. dm. Furthermore, we 
have the property tp*(t) = <p(t *), such that after the sub¬ 
stitution t = ix, Eq. m simplifies to 


p(s) 


pOO 

/ Re 

gixs e <p(ix) 

la 



dx. 


(14) 


Inserting Eq. (fT51) into this relation, we finally obtain 


e lT - 2(ir) 3 / 7 7 4/ 7(ir) + 3(ir) 2/7 75 /7 (ir) 


dr. 


(15) 
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Distributions p(s) obtained from numerical evaluation of 
Eq. (TTSl) are displayed in Fig. [5] as solid curves on top of 
the numerical results of Sec. EH The agreement is excel¬ 
lent for all the values of nz 3 . Furthermore, we notice that 
Eq. (1151) confirms the conclusion drawn from the simu¬ 
lations in Sec. IIII Dl p(s) depends on the parameter nz 3 
only, being completely independent of the microscopic 
details of the probe sphere and of the scatterers. 

D. Asymptotics 

Although the distribution p(s), Eq. (fl5l) . has no 
evident analytic expression, several simple asymptotic 
limits can be readily examined. 



FIG. 8. (color online) Distribution p(s) for nz 3 = 0.005 
(solid red curve), together with the approximate form (11911 at 
moderate s (dashed curve). 


Large-s limit We describe the limit of large s by expand¬ 
ing the term inside the square brackets in Eq. m up 
to second order in r < 1. Since essentially the val¬ 
ues of r such that rsnz 3 < 1 contribute to the inte¬ 
gral, this expansion is a good approximation provided 
1 /(snz 3 ) < 1. It results in a Gaussian integral which is 
straightforwardly performed to give: 


. . V33 nz 3 

V (s) ^ ^— exp 


33tt 3 2 


(16) 


Making use of Eq. m, we rewrite Eq. m as 


1 

, _ exp 

V 2ir5U 2 


(U-U) 2 ' 
2 SU 2 


(17) 


At large s, p(s) is thus simply a (normalized) Gaus¬ 
sian distribution of mean U and variance SU 2 . The 
expansion used to derive Eq. mi being valid as long as 
s ( nz 3 ) _1 , the Gaussian shape is a very good approx¬ 
imation of the whole distribution when nz 3 1, which 
is a direct consequence of the weakness of fluctuations 
at large distances. We expect this Gaussian distribution 
to be universal at large distances, as a consequence of 
the central-limit theorem (many scatterers contribute 
to U(z ) when nz 3 1), regardless the nature of the 
disordered medium. This conclusion is supported by 
similar predictions previously made in the context of the 
thermal Casimir effect in one-dimensional disordered 
media [25], [26|, as well as in recent studies of the CP 
interaction involving quasi two-dimensional disordered 
metals [27]. For comparison, we show Eq. m in 
the lower-right panel of Fig. [5] as a dashed curve, for 
nz 3 = 40. Note that in the chosen geometry this limit is 
difficult to reach in the numerical simulations, because 
it requires the generation of a significant number of 
scattering spheres in order to satisfy the condition 
nz 3 1, while maintaining z < A/2 to avoid finite-size 
effects. When nz 3 -C 1, Eq. (1171) still holds but only in 
the very far tail of the distribution, s 3> ( nz 3 ) _1 3> 1. 
The physical reason for which we recover a Gaussian 
tail even for small values of nz 3 is the following. Large 


values of s correspond to particular disorder realizations 
T> for which the Casimir potential is very large, i.e. for 
which the density of scatterers nx> underneath the probe 
atom is very high. Thus, for these specific disorder 
realizations the condition nx>~ 3 7S> 1 is effectively fulfilled. 


Moderate values of s To describe the limit s <C (nz 3 ) 1 , 
we expand the term inside the square brackets in Eq. 
m for large r: 


P (s) ~ -yj-j-Re 


dr exp 


27rnz J 


—isr — 2T(4/7) 


x (ir) 3/7 + 3T(5/7) (ir) 2/7 ] } . 


(18) 


We then perform a Wick rotation r = ix , and express 
both exponential terms exp(aa; 3 / 7 ) and exp(/3a; 2 / 7 ) as 
power series. This gives 


p(«) = ~lo Im 


. 7TUZ 

ax ex p (- sx 


n\m\ 


27TUZ 3 


E 

n,m=0 


T(4/7) (-*) 3/7 [7rn2 3 r(5/7) (-x) 2/7 


If we additionally assume s nz 3 , the terms ( m,n) = 
(0,1) and (m, n) = (1,0) give the leading contribution to 
p(s) (the term m = n = 0 is purely real and does not con¬ 
tribute). Keeping only these two terms and computing 
the remaining integral, we find 


p(s) 


2 1 n 3/7 O nz 3 ) 4 / 7 |" 1 f 7 xnz 3 s\ V ' 

7 s 10 / 7 V 10 ) 


(19) 


Equation (fTHl) holds for nz 3 <C s <C (nz 3 )- 1 . It is shown 
in Fig. [|as a dashed curve for nz 3 = 0.005, together 
with the exact distribution calculated from Eq. m 
(solid red curve). 

Small-s limit We finally consider the low-potential tail 
s < nz 3 of p(s). In order to find an asymptotic expansion 
in that limit, we come back to Eq. (USD and apply the 
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change of variables r = ( y/s ) 7 . This yields 




y 6 dy exp 


f 2nnz 3 
\ 6s 6 
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ly 


—2r(4/7)i 3/7 s 3 y 3 + 3T(5/7)i 2/7 sV } ■ (20) 

In the limit of very small s, we are thus led to evaluate 


I = Re dyy 6 exp [A f(y )\, (21) 

Jo 

where A = nnz 3 /(3s 6 ) 1 and f(y) = (3/10 )iy 7 — 

2r(4/7)i 3 / 7 s 3 y 3 + 3r(5/7)i 2 / 7 s 4 y 2 . Equation (l2ll) nat¬ 
urally calls for the method of steepest descent. There 
are five saddle points, solutions of f'(y ) = 0. Only one 
of them, denoted by y$p, turns out to give a nonzero 
contribution to p(s): 


V SP 


60T(4/7)s 

21 


3"I V 4 


e -i7r/14 


( 22 ) 


Using Cauchy’s theorem, we then deform the path of inte¬ 
gration to a path coinciding with the path of steepest de¬ 
scent in the vicinity of ysp ■ This is achieved by expanding 
f(y) up to second order around y = ygp and performing 
the change of variables x = {y — ysp) exp(—3i7r/14), leav¬ 
ing us with a Gaussian integral whose evaluation leads to: 


p(s) 




(23) 


with prefactors a = (20/7) 3 / 8 T(4/7) 7 / 8 and /3 = 
16 v / 27t( 5/7) 3 / 4 T(4/7) 7 / 4 /21. The asymptotic form (l23l) 
is completely analogous to the so-called Lifshitz tail that 
describes the band edge of the density of states of disor¬ 
dered conductors in solid-state physics |46|]. Physically, 
it can be understood from the following qualitative ar¬ 
gument. Low values of s are achieved for rare disorder 
realizations where the probe sphere stands above a large 
region free of scatterers, as illustrated in the right scheme 
of Fig. 0 Since the distribution of scatterers is Poisso- 
nian (Sec. IIVI) . the distance between the events of this 
Poisson process follows an exponential distribution. Con¬ 
sequently, the probability to find a large region of size R 
free of scatterers is oc exp(— cnR 3 ), where c is a numeri¬ 
cal constant. In such configuration, the Casimir potential 
felt by the atom can be estimated as (see Fig. 0: 


U = —2TrnT^ 


r(r — z)dr 


IVz' 2 +W 


oc 


—nT? 

R 4 


(24) 


to leading order in z/R <C 1. On the other hand, we have 
seen in Sec. IIV Bl that the average Casimir potential is 


U oc — nTi/z 4 . Therefore, for the rare disorder realiza¬ 
tion displayed in Fig. 0 we have U oc ( z/R) 4 U , such 
that 


p(s) ~ exp(— cnR 3 ) = exp 


U 


A/4' 


—cn 


' c / 1 / 4 


= exp 


s 3/4 


(25) 


which is nothing but the asymptotic form 


V. CONCLUSION 

In this paper, we have developed a statistical descrip¬ 
tion of Casimir-Polder potentials from both a numeri¬ 
cal and an analytical perspective. This approach is well 
suited for describing the Casimir interaction between 
a simple dielectric object and a strongly heterogenous 
medium made of a large number of independent consti¬ 
tuting elements. It can be readily extended to other ge¬ 
ometries and to heterogeneous media characterized by a 
more complex statistics involving, for instance, a non- 
uniform or polydisperse distribution of scatterers. 

As a first extension of our work, it would be interest¬ 
ing to investigate deviations to the dilute limit na 3 -C 1 
where the scatterers can no longer be systematically con¬ 
sidered independent. We expect these deviations to pri¬ 
marily affect the far tails of the distribution p(s). A 
second open question concerns the change in the statis¬ 
tics of the Casimir-Polder potential when the host (ho¬ 
mogeneous) medium has a dielectric constant differing 
from unity. This problem is more difficult to treat since 
it now involves a surface, which implies multiple reflec¬ 
tions inside and outside the medium. The presence of a 
host medium may therefore strongly affect the distribu¬ 
tion p(s) [in the obvious limit where the dielectric con¬ 
stant of the host medium goes to infinity, one recovers a 
perfectly reflecting interface and p(s) should tend to the 
Dirac function <5(s — 1)]. 

In practice, the distribution p(s ) could be experimen¬ 
tally accessed either by moving the sphere over a static 
disordered medium to record different disorder distribu¬ 
tions, or by taking advantage of a Brownian motion of 
the scatterers if the measurement process is fast enough. 
Indeed, in that case different disorder realizations can be 
obtained by detecting the Casimir force and then letting 
the scatterers move before carrying out the next mea¬ 
surement. If the measurement process is slow, the effect 
of the motion of the scatterers is to average the Casimir 
potential, giving him a value well approximated by Eq. 
0 since Doppler shifts have a negligible effect at thermal 
velocities I47l-l49l . 
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